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We point out some limits of the perturbative renormalization group used in statistical mechanics 
both at and out of equilibrium. We argue that the non perturbative renormalization group formalism 
is a promising candidate to overcome some of them. We present some results recently obtained in 
the literature that substantiate our claims. We finally list some open issues for which this formalism 
could be useful and also review some of its drawbacks. 
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I. INTRODUCTION 



Field theory together with renormalization group (RG) 
have provided very powerful means to investigate contin- 
uous phase transitions in equilibrium and non equilib- 
rium statistical physics. The need for field theoretical 
techniques is now well understood: when spatial (and/or 
time) fluctuations are large in a iV-body system, the 
mean field approach is inadequate and it becomes neces- 
sary to keep track of the space (and/or time) dependence 
of the order parameter. In the continuum, statistical fluc- 
tuations are summed over through a functional integral 
so that theories to be dealt with are (euclidean) field the- 
ories. The difficulty with these theories is that close to 
a continuous phase transition, fluctuations on all length 
(and/or time) scales contribute and render the perturba- 
tive approach problematic. The origin of the problem is 
well known and is twofold. First, at each order of the 
perturbation expansion (of, say, a correlation function) 
all fluctuations of all wave lengths (and/or frequencies) 
contributing to this order are summed over (Feynman 
graphs). Since they all contribute algebraically, the re- 
sulting sums (integrals) diverge at very short (ultraviolet) 
or very long (infrared) wavelengths. Getting rid of these 
ultraviolet divergences is the subject of perturbative 
renormalization and is now well under control^, 0, 0- 
Second, even after renormalization, the perturbation se- 
ries are in general non convergent and, as such, difficult 
to use to obtain reliable quantitative results. This is, for 
instance, the case for the perturbation expansion of the 
/3-functions of the RG flow and of the critical exponents. 
This is all the more severe that the dimensions of inter- 
est are far below the critical dimension since, then, the 
field theories are in their strong coupling regime. Fortu- 
nately, in some cases, the perturbation series turn out to 
be Borel summable so that efficient resummation tech- 
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niques (Pade-Borel, conformal mappings, etc) allow to 
compute physical quantities reliably and with high ac- 
curacy. This is the case for the 0(N) models in three 
dimensions for which five (e-expansion) and six (fixed di- 
mension expansion) loops have been computed and where 
all available methods lead to consistent and accurate re- 
sults for critical exponents^. Some of the most cele- 
brated equilibrium systems for which resummation tech- 
niques also work well are ferromagnetic systems with cu- 
bic anisotropy or with quenched disorder [a, la, 13]. How- 
ever, and contrary to the common belief, this ideal pic- 
ture turns out to be the exception rather than the rule. 
In most cases, the situation is more complicated for at 
least two reasons. First, in many cases, only the two first 
orders of perturbation series are known and are insuffi- 
cient to perform any resummation. This is the generic 
case in out of equilibrium statistical systems and for equi- 
librium ones involving fermions or gauge fields (as, for in- 
stance, the electromagnetic field). Second, even in cases 
where five or six loop series are known, resummations 
do not always lead to converged results. This is the 
case, for instance, for frustrated spin systems in three 
dimensions 0, 0, 0, HTl fl^L [l3| . Again, contrary to the 
common belief, the problem is not just to compute accu- 
rately critical exponents but to determine the qualitative 
behavior of the system at the transition: either first or 
second order phase transition, belonging or not to a given 
universality class, etc. 

Of course, apart from these difficulties, more serious 
ones — not only related to strong coupling behaviors — 
can exist and invalidate perturbation expansions. Many 
examples exist in different physical contexts. The low 
temperature expansion of the ferromagnetic XY (0(2)) 
model in two dimensions is identical at all orders to that 
of a free theory. This expansion is invalidated at finite 
temperature by the existence and liberation of vortices 
(Kosterlitz-Thouless transition) Q • Quantum Chromo- 
dynamics (QCD) is unable at any finite order of per- 
turbation to describe quark confinement. The pertur- 
bation expansion of the /3-function of the Kardar-Parisi- 
Zhang (KPZ) equation describing the growth of inter- 
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faces is known at all orders (it is given by the one-loop 
result) but is unable to describe the scaling behaviors 
in the rough phase above two dimensions [lH ITEj . The 
e = 6 — d expansion of the random field O(N) model 
(RFO(N)M) is identical at all orders to the e = 4 - d ex- 
pansion of the pure model (dimensional reduction) fifty . 
However, it has been rigorously proven that the RFIM in 
d = 3 and the pure Ising model in d = 1 display different 
critical behaviors [Til ITq . This implies that dimensional 
reduction must break down, at least below a dimension 
between six and three. 

In the two first examples, the problem is that the 
low energy excitations relevant for the description of 
the physics at large scale are qualitatively different from 
those at small scale: they are vortices for the XY model 
and bound states (actually confined states) like hadrons 
for QCD. Since, on the one hand, these excitations dom- 
inate the low energy physics and, on the other hand, 
they lead to non analytic contributions that are missed 
by perturbation expansions, the latters are unable to de- 
scribe the large scale physics of these systems. The pres- 
ence of non analyticities is also the reason for the break- 
down of dimensional reduction in the RFO(N)Mp3. It is 
not known, although suggested in the literature, whether 
these non analyticities are also related to the presence 
of a bound state|2£j |2l|, and neither it is for the KPZ 
equation. 

All these difficulties suggest that it is necessary to go 
beyond perturbation theory, at least to be able to deal 
with strongly coupled systems without having to resort 
to perturbation expansions at high orders and to rcsum- 
mations. Unfortunately, apart from two dimensions or at 
large N, no systematic method for solving exactly a field 
theory is known. One has to go back to the summation 
over the fluctuations in the partition function and try to 
organize it differently. This is what Wilson did with the 
momentum shell integration of rapid modes and the con- 
struction of effective theories for slow modes |22|. His idea 
was to integrate out fluctuations scale by scale and not 
order by order in a series expansion. Here, an important 
remark is of order. If it were possible to perform exactly 
this integration between two different scales, it would be 
possible to iterate this procedure until all fluctuations are 
integrated out. This would amount to solving exactly the 
problem and, of course, this is impossible in general. Ac- 
tually, there are only very few examples where RG tech- 
niques have enabled to solve exactly a model that was 
not solved (in a simpler way) by another method. In a 
sense, Wilson's procedure is more subtle. Of course, the 
momentum shell integration can be performed exactly in 
a formal way that is by expressing the RG flow of a quan- 
tity — e.g. an effective hamiltonian — in terms of itself. 
This leads to the formulation of an exact RG equation 
(that we derive in the following and which is a functional 
differential equation) [2^, [2j| . But this equation cannot 
be solved in general. The important point in Wilson's 
method is that the very idea behind it suggests a natural 
way to approximate the integration over the rapid modes 



and thus a natural truncation of the exact RG equation. 
Thus, RG techniques arc almost intrinsically linked with 
approximation methods: they are not very interesting 
when exact solutions are known. 

The aim of this article is to show how Wilson's RG 
ideas can be concretely implemented both at and out of 
equilibrium. In particular, we shall show how Wilson's 
idea of a flow of effective hamiltonians is conveniently re- 
placed by that of a flow of Gibbs free energies and we 
shall generalize it to non equilibrium statistical mechan- 
ics where the free energy no longer exists. One of the cru- 
cial part of this article will, of course, be devoted to the 
most developed truncations of this flow. In particular, 
we shall show that they allow both to accurately deter- 
mine universal quantities, such as critical exponents even 
when the theories are strongly coupled, and non univer- 
sal quantities, such as phase diagrams. This last point 
is a crucial advantage of this method as compared with 
perturbation theory. 



II. THE NON PERTURB ATIVE 
RENORMALIZATION GROUP (NPRG) 
FORMALISM 

A. The equilibrium case 

Any implementation of Wilson's momentum shell in- 
tegration relies on a separation between rapid and slow 
modes. This separation can be achieved through a sharp 
cut-off in momentum space but this procedure leads to 
singularities when approximations are performed. It is 
preferable to perform a smooth separation by modify- 
ing the Boltzmann weights in such a way that the slow 
modes effectively decouple from the model while the 
rapid ones are unaffected. This is conveniently achieved 
by adding to the original partition function (for equilib- 
rium systems) a "momentum-dependent mass term" giv- 
ing a large "mass" to slow modes and leaving unchanged 
the rapid ones. The slow modes are thus weakly corre- 
lated, they do no longer propagate and do not contribute 
to the long distance physics of the (modified) model. For 
a theory with one scalar field <f> and hamiltonian H[cj)], 
we hence define a scale dependent family of partition 
functions^ |H |H IH |H S3 : 

Z k [B] e-^-^M+Z B -\ (1) 

where B is an external source, e.g. a magnetic field, k 
a momentum scale and A.Hk[(f>] is the scale dependent 
mass term responsible for the decoupling of slow modes 
in Z k : 

AH fe \4>\ = \J d d x R k (x - y)cp{x)cp{y) (2) 
= \J (0ffl*(9WffM-3). (3) 
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R k is the cut-off function that must verify: 



q 2 » fc 2 : R k (q 2 ) 
q 2 « fc 2 : R k (q 2 ) 



> no mass for rapid modes (4) 
k 2 large mass for slow modes, (5) 



in order to achieve the decoupling. A typical and very 
simple function satisfying these constraints is|30j| 



Rk(q) = (fc 2 - q 2 )9(k 2 - q 2 ) 



(6) 



which is a parabola with negative concavity on [0, fc] and 
vanishes elsewhere. These definitions yield that: 

1) when k = 0, Rk=o(q) = identically and thus 
Z k= o = Z, the original partition function of the model 
under study. The latter is therefore recovered in the limit 
fc — * 0, that is when all fluctuations have been integrated 
out. 

2) when fc = A ~ inverse lattice spacing (or, more 
generally, the microscopic scale), Rk{q) is very large for 
all modes so that all fluctuations are frozen. Thus 
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MF 



log Z k= 

(7) 

where M MF is the configuration with highest Boltzmann 
weight, that is M MF is the mean field magnetization of 
the model with hamiltonian H + AH\. 

Then, if one defines the family of functionals r/-[M] 
given by the Legendre transform of the free energy 
Wfc[_Bl = log Z k [B] (up to the last term proportional to 

R k )m- 

T k [M] + W k [B] = JBM-^J R k (q 2 )M q M_ q (8) 



with 



M(x) 



SW k 
5B(x) ' 



(9) 



it follows that: 

1) when fc = 0, T k —Q = T is the Gibbs free energy 
of the system since it becomes in this limit the standard 
Legendre transform of W (R k= o is identically vanishing) ; 

2) when k = A, F k=A [M] ~ H[M] from Eqs.0 and 
©. 

Thus, r^[M] smoothly interpolates between the micro- 
scopic hamiltonian and the Gibbs free energy when fc is 
decreased from A to 0, that is, when more and more flu- 
cuations are integrated out[2!|. The interpretation of Tk 
and M is simple. When B(x) is kept fixed, M as defined 
by Eq.@ is fc-dependent. It is the precursor at scale fc of 
the physical magnetization and corresponds to the mag- 
netization of a system composed of blocks of spins of size 
fc -1 , that is of a system where only fluctuations inside 
blocks of size fc -1 have been summed over. rj,[M] is the 
Gibbs free energy for the configuration with given mag- 
netization M of this system. Thus, contrary to Wilson's 
effective hamiltonian at scale fc, which is the hamiltonian 
for the slow modes that have not yet been integrated out 



in Z, T k [M] is the free energy for the rapid modes that 
have already been integrated out. It enables, in principle, 
to compute correlation functions of the rapid modes, con- 
trary to Wilson's effective hamiltonians. I^fM] is called 
the effective average action (rfc = o[Af] is called the effec- 
tive action in the context of quantum field theory; it is the 
generating functional of the 1-Particle-Irreducible (1PI) 
correlation functions) . 

The exact evolution equation of T/JM] with fc is ob- 
tained by differentiating Eq.Q with respect to fc and by 
eliminating W fc for T k [M] thanks to Eq.©. It writes [29|: 



1 



d s T k = - I R k (q 2 )(r k 2 \q,-q) + R k (q 2 )) 



where s = log(fc/A), R k = d s R k , 



T{ n \q u ...,q n7 M) = 



S n T k [M] 



SM qi . . . 6M qn 



(10) 



(11) 



and (Tjp(q, —q) + R k (q 2 )) 1 is the inverse, in the oper- 
ator sense, of T k 2 \q,—q,M) + R k (q 2 ). Eq. lJTUj l is exact 
and, as such, contains all perturbative and non perturba- 
tive physics. It is a functional (since T k ' depends on M), 
partial differential equation that requires truncations to 
be solved. The very idea of incomplete integration and 
effective theories is to preserve as much as possible the 
long distance physics, approximating the short distance 
one. Thus, it is natural to expand r^[JVf] in powers of 
derivatives of M. This is the so-called derivative expan- 
sion. In practice, this amounts to proposing an ansatz for 
[M] which is polynomial in the derivatives of the order 
parameter. In the case of a one component scalar field 
theory with IL^ symmetry, this consists in considering for 
instancelUlii: 
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d d x U k {p) + -Z k (p)(VM) 2 + 0(V 4 ) 



T k [M]-- 

(12) 

where p = M 2 /2, U k and Z k (that should not be confused 
with the partition function Z k ) are functions of M 2 and 
fc. For uniform magnetization, T^fM] reduces to U k (up 
to a volume factor) which is called the effective potential. 
The term proportional to Z k (p) is the first correction to 
U k taking into account the momentum dependence of 
the correlation functions T, n \qi, . . . , q n ). Z k (p) enables 
to compute the anomalous dimension of the field M. Of 
course, the ansatz l)12|l cannot be an exact solution of 
r/-[M] and thus U k (p) and Z k (p) must be defined in such 
a way that their flows do not get any contribution from 
the 0(V 2 ™) terms with n > 1 that are supposed to be 
neglected in the ansatz. It is convenient to define them 
as: 



U k {p) 



r(p) 
v 

i d 



^ = ^r<f>(«,- ff )(p)| 



(13) 
(14) 
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where V is the volume of the system and p is a uniform, 
that is x-independent, field configuration. Inserting these 
definitions in the RG flow equation (fTT)f) leads to the RG 
equation for Uk{p) and Zk(p). These equations are par- 
tial differential equations. It is interesting to notice that 
it is possible to further truncate while preserv- 

ing qualitatively and semi-quantitatively the properties 
of the method. This additional truncation consists in 
field expanding Z k {p) and/or the potential Uk{p), that 
isHEa: 

1) Uk{p) is truncated by keeping only a polynomial in 
p of low degree. It is convenient (and more accurate) in 
this case to expand Uk(p) around its running minimum 
Kk- The simplest non trivial field truncation is: Uk(p) = 

\ k /2{p- K k f +0{{ P - K k f). 

2) Zk(p) is truncated by keeping only a field- 
independent coefficient: Zf.(p) —tZ^ = Z k {p = k^). 
This is the simplest truncation that allows to compute 
a non trivial anomalous dimension. Let us notice that it 
has for long been believed that the anomalous dimension 
was essentially arbitrary within the derivative expansion 
and that it was therefore useless to include a %(p) / 1. 
Taking Zk{p) — I — and thus r\ — — is called the Local 
Potential Approximation (LPA) and was first used (in a 
slightly different way) by Wegner and Houghton back in 
the seventies |3l|. 

If we combine truncations 1 and 2, the RG flow be- 
comes a set of three simple coupled differential equations 
for Kfc, Afc and Zfc[25|]. A very important point with this 
truncation when it is implemented on the O(N) mod- 
els is that it allows to recover i) the one-loop result ob- 
tained perturbatively in d = 4 — e, ii) the N — > oo result 
in any dimension and, for N > 3, iii) the one-loop re- 
sult obtained in d = 2 + e from the non linear sigma 
model0,|2jj. This is very different from the perturba- 
tive results that are very difficult to extrapolate beyond 
their domain of applicability. Let us also notice that for 
N = 2 and d — 2, this truncation enables to retrieve most 
of the qualitative as well as some quantitative features of 
the Kosterlitz-Thouless transition|32| . 

There are two origins for the errors made with an 
ansatz such as Ea. ljT2"|l supplemented or not by a field 
expansion of Uk{p) and/or Z k {p). First, the exact RG 
equations for Uu and Z k involve the neglected terms of 
the derivative expansion of order V 4 , V 6 , .... Moreover, 
if a field expansion of Uk(p) is performed as in trunca- 
tion 1) above, the RG equations of n k and involve 
the coupling constants of the next terms of the expan- 
sion, that is the coefficients of (p — Kfc) n with n > 2. 
Second, the choice of a cut-off function R k that, in prin- 
ciple, has no influence since Rk(q 2 ) vanishes identically 
in the limit k — > 0, does matter once truncations are per- 
formed. Many studies have been devoted to finding an 
optimal choice |33l 133. l35l |3(| |37| None of them gives a 
complete solution to this problem. 

Let us now work out the RG equations obtained with 
a standard truncation as l|12|) for the effective potential 
and anomalous dimension. Let us first recall that since at 



a second order phase transition the correlation length is 
infinite, the system becomes scale invariant at the tran- 
sition point (at sufficiently large distance compared with 
the lattice spacing). Thus, in units of the (inverse) run- 
ning block spin size k, it becomes invariant under the RG 
transformations (for t « A). It is therefore convenient to 
work with dimensionless quantities: in terms of them the 
transition corresponds to a fixed point of the RG flow. 
Actually, to find out a fixed point, it is also necessary to 
work with a renormalized magnetization. Indeed, at the 
transition point, the magnitude of the running (sponta- 
neous) magnetization M never stops flowing towards as 
k — ► 0, which a priori prevents from finding a fixed point. 
It is thus necessary to work with a (dimensionless) renor- 
malized magnetization whose magnitude tends to a finite 
value as k — > 0. We thus define 



P - 

Uk{p) - 



z k k 2 - d P 

k- d U k (p). 



(15) 
(16) 



With truncation 2) above and with the regulator of 
Eq.lJBJ, the RG equations for the potential Uk(p) and for 
Zk become very simple: 



d s u k {p) 



-- ~duk{p) 

(l- 
d \ d 



(d-2- 
Vk 



Vk)pu' k (p) 



1 



Vk 



-d s log Z k = 



+ 2 
I6v d 



!+<((>) 



2K(P) 



d 1 + 2XkKk 



(17) 



(18) 



where Vd — l/(2 d+1 7r rf / 2 r(d/2)), primes denote deriva- 
tives with respect to p, Kk is the minimum of u k (/?) and 
Afe = u^(Kfc). rjk is the running anomalous dimension 
which goes to a fixed point value as k — > — the physi- 
cal anomalous dimension — if the potential itself reaches 
a fixed point solution. 

This theoretical framework — or those deriving from 
more sophisticated truncations — have turned out to be 
very powerful. We show in the following some results 
that have been obtained for the Ising and 0{N) mod- 
els by ourselves and by other groups. However, before 
this, let us show how to adapt this formalism for out of 
equilibrium statistical systems. 



B. The out of equilibrium case 

For out of equilibrium systems the very notion of free 
energy does no longer exist in general. However, it is still 
possible to define generating functionals of correlation 
functions, analogous to the partition function Z. More- 
over, even if it is no longer possible to speak of Boltzmann 
weights and hamiltonians, the generating functionals Z 
are still often given by functional integrals of weights 
that, formally, enable to define "hamiltonians", called ac- 
tions in this context. The Legendre transform T of log Z 
is then analogous to the Gibbs free energy and is the 
generator of the 1PI correlation functions. There are two 



5 



main means to build field theories in out of equilibrium 
statistical mechanics: either from a Langevin formulation 
of the problem or from a more microscopic approach — 
if it exists — a la Doi-Pelitij3^,|3!|, that is from a master 
equation. 

Even if approaches at and out of equilibrium arc close 
together, some physical differences are crucial, such as 
the irreversibility of the dynamics, the violation of the 
fluctuation-dissipation theorem, etc. Technically, these 
differences have a translation in the formalism: i) the 
field theories are written not only in terms of the "phys- 
ical" field(s) but also in terms of the response field(s), 
ii) temporal fluctuations are taken into account by a sin- 
gle time derivative (contrary to the Laplacian for spatial 
fluctuations), iii) the actions are not "hermitic", that is, 
the physical and response fields do not necessarily play 
symmetric roles 40]. However, all these differences do not 
prevent from constructing an analog of the effective aver- 
age action method for out of equilibrium svstems |4lll42| . 
The non trivial point is to choose a R k function. Should 
we take it space and time dependent? acting only on 
the physical field or on both the physical and response 
fields? Since the free field theory corresponds to the fol- 
lowing theory: 

Z = J D(f>D4>e- So (19) 
with (20) 
So = J d d xdt 4>(x, t) (d t - DV 2 ) <j>{x, t) (21) 

it is natural to define the Rk term by: 

ASk = J d d xd d ydtdt' $(x,t)R k (x - y,t - t')<j)(y,t'). 

(22) 

We have also chosen to take it time-independent so that 
the time integral in the RG equation can be trivially 
performed|4l|. It is not known whether this corresponds 
to a "good" choice from the accuracy point of view. Once 
these choices have been made, it is possible for a non triv- 
ial model with action S to define: 

Z k [J,J]=J 2ty2ty e -S[^]-AS* ! [<M+/ *4>+J 34 _ ^ 
The effective average action is defined as in Eq.© by: 

r fe [v, $ + w k [J, J] = / JV' + / H (24) 

R k (q 2 )ip q , u rp- q ,- u (25) 

J q 

with 



#M) = ^ and tf( M )= ^ 



5J(x,t) 



5J{x,t)' 



(26) 



to if> and i> and the 2x2 matrix Rk, the off-diagonal ele- 
ments of which are Rk (q 2 ) ■ Then, the out of equilibrium 
NPRG equation writes [4l| : 



~sR k (f k 2) (q,Gj;-q,-u) + R k (q') 
A typical ansatz for a reaction-diffusion problem is: 

r k [ip,i>}= [ d d xdt {u k {ip^) + z k ^{dt-D k v 2 )^} 



(27) 



Within this formalism some problems of branching 
and annihilating random walks (BARW) have been 
addressed^! that we briefly sketch further on. 



III. SOME RESULTS OBTAINED WITH THE 
NPRG METHOD 

Among others, a crucial problem as for the NPRG 
method combined with the derivative expansion, is to es- 
timate the quantitative accuracy of the results obtained. 
More precisely, it would be highly desirable to be able to 
estimate error bars. As usual, this is a very difficult task 
when no analytic result is known on the behavior of the 
series encountered in the method. As already mentioned, 
this question is closely related to the choice of the cut-off 
function Rk- 



order 


v rj 




0.6506 


d 2 


0.6281 0.044 


d 4 


0.632 0.033 


7-loops 


0.6304(13) 0.0335(25) 



A (2) 

II is convenient to define both the 2x2 matrix T y k [tjj,tp] 
of second functional derivatives of [ip, ip] with respect 



TABLE I: Critical exponents of the three dimensional Ising 
model. d°, d 2 and d 4 correspond to the order of the trun- 
cation of the derivative expansion (NPRG method) |36j . For 
completeness, we have recalled in the last line, the results 
obtained perturbatively|4|. 



We have studied the critical exponents of the Ising 
model in three dimensions as a testing ground of the con- 
vergence of the derivative expansion at order V° (LPA), 
V 2 and V 4 36] . The rule of thumb adopted consists in 
evaluating the error bar through the evolution of the val- 
ues of the exponents with the order of the truncation. In 
the different cases studied, we have found that the dis- 
tance between the best known results and ours decreases 
when increasing the order of the truncation of the deriva- 
tive expansion. Thus, it seems reasonable to believe that 
the error made at a given order of the derivative expan- 
sion is given by the distance between the results of the 
last two orders. More precisely, at each order of the trun- 
cation, we have studied the dependence of the critical ex- 
ponents on the choice of function Rk- Each time, we have 
tried to find the best R k by requiring the exponents to 
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depend as weakly as possible on the choice of Rk (since 
they verify this property in the exact theory). We have 
therefore considered a one-parameter family of functions 
Rk and have computed the critical exponents for all el- 
ements of this family. The function Rk that is selected 
as optimal is that for which the exponents are stationary 
with respect to a change of Rk (this is an implementa- 
tion of the Principle of Minimal Sensitivity) [35| ■ At each 
order of the derivative expansion, this leads to a set of 
optimal exponents. At order V , these optimal expo- 
nents are quite close to the best known values, see Table 
I. Thus, although only three orders of the derivative ex- 
pansion are known, it seems that, at least in the three di- 
mensional Ising case, the derivative expansion converges 
rapidly. It is important to notice that no resummation 
of any kind has been necessary to obtain these results. 

As a conclusion, one can see on Table I that the results 
obtained from the NPRG method are in a fairly good 
quantitative agreement with the best known results for 
the three dimensional Ising model. This is, of course, very 
encouraging as for the reliability of the method and the 
convergence of the derivative expansion. Let us moreover 
emphasize that the Ising model belongs to the very small 
class of models, the perturbation expansion of which in- 
volves only one coupling constant and is known to be 
Borel summableQ. In many cases, things are much more 
complicated. A famous example is given by the three 
dimensional frustrated spin systems, either the triangu- 
lar antiferromagnets or the helimagnets. Here, the re- 
sults obtained by means of the resummed perturbation 
series and those obtained from the NPRG are qualita- 
tively different d U [H EJ 13 El The former predicts 
a second order phase transition whereas the latter pre- 
dicts very weakly first order transitions. The debate is 
still open and it will be interesting to see which prediction 
is correct. 

Let us now turn to out of equilibrium systems and more 
precisely to branching and annihilating random walks. 
The problem is the following. On a d-dimensional lat- 
tice, particles are diffusing (with rate D) and can create 
offsprings through the reaction A — > A + A, with rate 
a. When they meet on the same site, they can annihi- 
late through A + A — > with rate A. The problem is 
to determine whether the system undergoes a continuous 
phase transition between an active and an inactive (ab- 
sorbing) phase and to which universality class it belongs. 
The field theory is obtained from the master equation 
and the action reads Ha: 



S = d d x dt {4> (d t - DV 2 



/ 2ctA(00 2 - 4> 2 <t>) 
\4> 2 <j) 2 }. (29) 



First of all, a mean field study predicts that the system 
is always in the active phase for a =/= 0. From a pertur- 
bative analysis performed at and below d = 2, Cardy 
and Tauber have obtained the phase diagram for small 
X/D and cr/D, showing unambigously that fluctuations 



are able to destabilize the active phase so that an ab- 
sorbing phase exists at finite a/D for d < 20. As 
for the universality class, this phase transition belongs 
to the Directed Percolation one. However, their calcu- 
lations also suggest that for all dimensions larger than 
two, the fluctuations are not strong enough to destabi- 
lize the mean field result and that, therefore, the system 
is always in the active phase. Using the NPRG method, 
we have re-examined the physics of this system. We have 
been able to compute the critical exponents in all dimen- 
sions and have shown that the phase diagram obtained 
perturbatively is wrong^. More precisely, above di- 
mension two, the system is indeed always in the active 
phase at small a/D and X/D as predicted by the pertur- 
bative analysis. But there is a threshold value of X/D 
above which an absorbing phase can exist. By numer- 
ically integrating the RG flow from dimensions 3 to 6 
we have determined these threshold values as well as the 
complete phase diagram 42J. In Figure ^ is displayed 
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FIG. 1: Phase diagrams of BARW A ^* 2 A, 2 A A in 
dimensions 1 to 6. Lines present NPRG results, rescaled as 
explained in the text. Symbols follow from numerical simula- 
tions. For each dimension, the active phase lies on the left of 
the transition line, the absorbing phase on the right. 

both the NPRG results and those obtained from numer- 
ical simulations. Only one free parameter has been used 
to match the numerical and analytical results. This pa- 
rameter takes into account the fact that both approaches 
differ at the microscopic scale A since one is defined on 
the lattice and the other in the continuum. However, it is 
known that both models should match after some RG it- 
erations since, at sufficiently large scale, the lattice model 
becomes homogeneous and isotropic and therefore can be 
described by a continuum theory. Thus, after a proper 
matching of the ultraviolet scales A of both models - 
requiring one rescaling parameter — the results can be 
compared quantitatively. The agreement is excellent and 
proves that, also in this example, the NPRG method is 
reliable. Let us emphasize that the phase diagram is a 
non universal feature of the model. Computing it requires 
to keep track all along the RG flow of both relevant and 
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irrelevant couplings since they all contribute to the non 
universal behavior of the model. While this is in general 
a formidable task in perturbation theory, that anyway 
is useless if these couplings are large, it just consists in 
the NPRG approach in integrating the flow equations for 
the potential with the bare potential of Eq. I|29(l as initial 
condition at scale A |43. E3. 14^ . 



IV. CONCLUSION AND SOME OPEN 
PROBLEMS 

Let us now come to a conclusion. The NPRG method 
seems to be a promising way to tackle the strong coupling 
regime of (at least some) field theories while avoiding 
problematic resummation methods used in perturbative 
renormalization. It also enables to compute non universal 
quantities such as phase diagrams because it can handle 
the functional nature of the effective potential with both 
its relevant and irrelevant couplings, that are most of the 
time infinitely manv|42l |45|. For the reaction-diffusion 
problem we have discussed in this article, the NPRG 
has allowed us to find that a phase transition exists in 
all dimensions while perturbative renormalization failed, 
probably even in a all-order analysis. This ex amp le (and 
it is not unique, see for instance [L^l4^l4alT7ll43 . l49| ) al- 
ready shows that this method leads to highly non trivial 
results. Now, an interesting question is: can we clas- 
sify problems for which perturbative theory could fail 
whereas this method could still be relevant and, recipro- 
cally, can we list its weak points? Let us try. 

We have already mentioned the problems encountered 
in the three dimensional frustrated systems for which the 
NPRG and the six-loop perturbative computations lead 
to different resultsfl^. It is interesting to understand 
why, since this problem probably belongs to a whole class 
of problems where perturbative theory is problematic and 
where the NPRG method could be fruitful (see H3)- 
The problem is the following. The Ginzburg - Landau 
- Wilson field theory (</> 4 -like) relevant for the descrip- 
tion of the critical physics of a system is determined by 
three properties of the system: i) the symmetry breaking 
scheme, G — > H, between the high (G) and the low (H) 
temperature phases, ii) the number N of components of 
the order parameter and iii) the space dimension d. The 
number of coupling constants in the hamiltonian depends 
on the symmetry group G and the lower this symme- 
try, the larger the number of invariant couplings in the 
hamiltonian. For all systems with a symmetry lower than 
O(N), at least two coupling constants are present in the 
hamiltonian. In these cases, the number of fixed points is 
larger than one, contrary to the O(N) case where there is 
only one. When N and/or d are varied, the fixed points 
move in the coupling constant space and it is generically 
observed that, at fixed d, there exists a critical value N c 
of N where two of them collapse. Two possible situa- 
tions can then occur: either both fixed points disappear 
when N is further varied (they become complex) or they 



exchange their stability. These two situations are, for 
instance, encountered in frustrated systems and in ferro- 
magnetic systems with cubic anisotropy. The difficulty is 
that as d is varied, N c changes and, of course, the prob- 
lem is to compute N c (d) since it is crucial to determine 
if the fixed points exist in d = 3 for N = 1,2,3 and, if 
they do, which one is stable. This is a highly non trivial 
task in perturbation theory that requires, in the known 
cases, resummations of the series obtained at least at 
four loops. Clearly, if the NPRG method leads to quan- 
titatively reliable results in these cases, this would avoid 
tedious perturbative calculations, that are anyway use- 
less if more than six loop results are necessary to obtain 
a converged value of N c (d = 3). These problems belong 
to the whole class of "strong coupling problems" of field 
theory. 

We have already mentioned another class of prob- 
lems where NPRG could be useful: those where per- 
turbation theory fails at all orders to describe the crit- 
ical physics. An important example is given by the 
RFO(N)M where dimensional reduction has been proven 
to be wrong [T7I Hsj . Precisely on this example, the 
NPRG approach has been used very recently [lQ]. The 
origin of the problem is that the effective potential Uk 
acquires non analytic contributions at finite k that are 
missed by perturbation theory. These non analyticitics 
have been successfully computed from the NPRG method 
since the functional nature of the effective potential can 
be dealt with. It will be very interesting to see if this 
method is able to obtain qualitatively and quantitatively 
correct results for the Ising case in three dimensions. 
Anyhow, the present results in this model already open 
the way to a renormalization group which is both func- 
tional and non perturbative. 

A third instance where NPRG could turn out to be 
useful is the study of phase transitions where topological 
defects play a role (especially in low dimensions). Let 
us recall that almost nothing is known in the physics of 
phase transitions induced by vortices in two dimensions 
apart from the Kosterlitz-Thouless transition for which 
the decoupling between spin waves and vortices allows 
a perturbative study of the transition. It could seem 
strange at first sight that, in situations where non trivial 
spatial configurations of the field are crucial, anything 
relevant can be learnt from a derivative expansion that 
keeps only the lowest order(s) in the momentum depen- 
dence of correlation functions. Let us however make two 
general remarks. First, the Kosterlitz-Thouless transi- 
tion has been studied by the NPRG method without in- 
troducing the vortices by hand d la Villain. Almost all 
the physics is quantitatively recovered apart from the fact 
that the line of fixed points is replaced by a line of quasi- 
fixed points where the /3-function of the temperature is 
very small although not strictly vanishing 48) . Second, 
spatially non trivial field configurations can play a great 
role in the summation over fluctuations, that is for the 
computation of the functional integral, while correlation 
functions, that are the results of this summation, can be 
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very smooth functions of the momenta. Thus, it is possi- 
ble that a derivative expansion succeeds in capturing the 
physics of topological defects. This would be extremely 
interesting since in many cases other than the 0(2) case, 
no Villain trick is available to write down a field theory 
that can be studied perturbatively. 

Let us now review some of the drawbacks of the NPRG 
method that would clearly require much effort to be over- 
come. First, the derivative expansion forbids to compute 
the momentum dependence of correlation functions. It 
just focuses on the flow of coupling constants with em- 
phasis on those of the effective potential. It is there- 
fore insufficient to compute cross sections in particle 
physics or structure factors. Second, it is not easy to 
point out the "small parameter" that underlies the va- 
lidity of the derivative expansion. We can postulate that 
the larger the anomalous dimension, the worse its con- 
vergence. However, this important assumption has yet 
not been supported by a systematic theoretical analysis. 
Third, the very way the cut-off function Rk(q 2 ) is intro- 
duced can be problematic. For instance, it breaks gauge 
symmetry and it is clearly one of the big challenge of the 
NPRG to allow gauge invariant computation jH^]. Fourth, 
the choice of an ansatz for Tk relies on the choice of its 
field content. In statistical mechanics, contrary to parti- 
cle physics, the microscopic degrees of freedom, at scale 
A, are known in general and the difficulty is to propose 
an ansatz in terms of the low energy degrees of freedom. 
As already emphasized, they can qualitatively differ from 
those at high energy — e.g. they are bound states — and 
it is non trivial to know whether the RG flow of, say, a 
potential will signal this change. The problem is that, 
although the physics can still be described in terms of 



the microscopic degrees of freedom, it can then turn out 
to be very complicated. For instance, while the physics 
is "local" in terms of the bound state field, it can become 
non local in terms of the microscopic degrees of freedom. 
Thus, a derivative expansion formulated in terms of the 
bound state field can be accurate whereas it is not in 
terms of the original fields [HJ. The choice of low energy 
degrees of freedom is therefore crucial for devicing an 
ansatz for and it is not clear up to now if there exists 
a systematic way for, at least, detecting a bad choice. 

NPRG is still in its infancy. For sure, some of the 
above mentioned problems will be solved in the future. 
Some others will remain out of reach. As usual, numer- 
ous physical problems will challenge us and will force us 
to better understand the NPRG method as well as its 
intrinsic limits. 
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